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^ ; Abstract 

A forefront area of research concerns the exploration of the properties of hadronic matter under 
extreme conditions of temperature and density, and the determination of the equation of state-the 
relation between pressure, temperature and density-of such matter. Experimentally, relativistic 
heavy-ion collision experiments enable physicists to cast a brief glance at hot and ultra-dense 
matter for times as little as about 10~ 22 seconds. Complementary to this, the matter that exists 
in the cores of neutron stars, observed as radio pulsars, X-ray pulsars, and magnetars, is at low 



2 ' temperatures but compressed permanently to ultra-high densities that may be more than an order 

of magnitude higher than the density of atomic nuclei. This makes pulsars superb astrophysical 



laboratories for medium and high-energy nuclear physics, as discussed in this paper. 



x. 

o3 1 Introduction 

Exploring the properties of hadronic matter under extreme conditions of temperature and density has 
become a forefront area of modern physics, both theoretically and experimentally. On the earth, heavy 
ion experiments enable physicists to cast a brief glance at such matter for times as little as about 
10 -22 seconds. On the other hand, it is well known that galaxies like our Milky Way contain up to 
10 8 to 10 9 neutron stars, which are observed as pulsars (rotating neutron stars). Such objects contain 
ultra-dense hadronic matter as a permanent component of matter in their centers. Radio telescopes, 
X-ray satellites-and soon the latest generation of gravitational-wave detectors-provide physicists and 
astronomers with an unprecedented wealth of high-quality data on such objects and thus serve as the 
observational windows on the inner workings of pulsars. This feature makes pulsars superb astrophysical 
laboratories for medium and high-energy nuclear physics (plus other fields of physics) [IlEllHllllGSElEllE]- 
Some of the key questions that can be addressed by studying the properties of pulsars are (see Fig. [Q: 

• What are the fundamental building blocks of cold ultra-dense matter? Specifically, does "exotic" 
matter exist in the cores of pulsars, such as boson condensates, color-superconducting quark 
matter, and multi-quark states? 
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Figure 1: The multifaceted connection between compressed baryonic matter matter and 
pulsar physics [Bj. 

• Are there pulsar observables that could signal the existence of exotica in their cores? 

• Are there rotationally-driven (accretion-driven) phase transitions in pulsars? 

• How does color-superconducting quark matter alter the properties of pulsars? 

• What is the true ground state of the strong interaction? Is it ordinary nuclear matter (i.e. atomic 
nuclei) or a color-neutral collection of up, down, and strange quarks (so-called strange quark 



• What are the distinguishing features of pulsars made of strange quark matter rather than ordinary 
hadronic matter? 

• What are the properties of matter subjected to ultra-high density radiation fields, ultra-high 
magnetic fields, ultra-high electric fields? How do such fields alter the properties of pulsars? 

• What are the key nuclear (heavy ion) reactions in the non-equilibrium crusts of accreting X-ray 
pulsars? 

• How strongly do pycnonuclear reactions in the crusts of accreting neutron stars alter the thermal 
evolution of such objects? 

• Do gravitational-radiation reaction driven instabilities limit the spins of pulsars? 

• What is the shell structure for very neutron rich nuclei in the crusts of pulsars? Do N=50 and 
N=82 remain magic numbers? 

Before we discuss some of these issues in more detail, we review some of the key properties of neutron 
stars. The first property concerns the masses of neutron stars, which range theoretically from around 
0.1 M Q (where M & = 2 x 10 33 g is the mass of the sun) to about 3 M & . Matter in the centers of neutron 
stars possess densities ranging from a few times n to an order of magnitude higher. Here n = 0.15 
nucleons/fm 3 denotes the baryon number density of normal nuclear matter, which corresponds to a 
mass density of 2.5 x 10 14 g/cm 3 . The number of baryons forming a neutron star is of the order 
of 10 57 . Rotating neutron stars are called pulsars. Three distinct classes of pulsars are currently 
known. These are (1) rotation-powered pulsars, where the loss of rotational energy of the star powers 
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Figure 2: Models for the EoS (pressure versus energy density) of neutron star matter [5J. 
The notation is as follows: RMF=relativistic mean-field model; DD-RBHF=density de- 
pendent relativistic Brueckner-Hartree-Fock model; n=neutrons; p=protons; H=hyperons, 
K=K~[u,s] meson condensate; Q=u,d,s quarks; H-matter=H-dibaryon condensate. 



the emitted electromagnetic radiation, (2) accretion-powered (X-ray) pulsars, where the gravitational 
potential energy of the matter accreted from a low-mass companion is the energy source, and (3) 
magnetars (e.g, SGR 1806-20), where the decay of a ultra-strong magnetic field powers the radiation. 
The fastest, very recently discovered neutron star, PSR J1748-2446ad, rotates at a period of 1.39 ms 
(which corresponds to 719 Hz) [9J, followed by PSR B1937+21 [TO] and PSR B1957+20 p] whose 
rotational periods are 1.58 ms (633 Hz) and 1.61 ms (621 Hz), respectively. PSR B1957+20 is moving 
through the galaxy at a speed of almost a million kilometers per hour. Due to this remarkable motion 
a bow shock wave is visible to optical telescopes. Other significant pulsars are Cen X-3 (first X-ray 
pulsar), SAX J1808.4-3658 (first accreting millisecond X-ray pulsar), and PSR J0737-3039A&B (first 
double pulsar binary system, which will permit a strong- field test of General Relativity). When neutron 
stars are formed they have interior temperatures of the order of 10 11 K (around 10 MeV). They cool by 
neutrino emission processes to interior temperatures of ~ 10 10 K within a few days. Throughout most 
of the active life of neutron stars in pulsars and X-ray sources, the interior temperature is ~ 10 7 — 10 9 
K. Surface temperatures are an order of magnitude or more smaller. Measurements of the surface 
temperature of the Crab pulsar (B0531+21), for instance, have led to T < 2 x 10 6 K.) Rotating 
magnetized neutron stars in pulsars are surrounded by a plasma, the so-called magnetosphere, in which 
via plasma processes the electromagnetic radiation from pulsars is generated. The magnetic fields near 
the surface of neutron stars in compact X-ray sources play an important role in channeling the accreting 
matter onto the neutron star surface. To date over 1600 pulsars are known. This figure is expected to 
increase dramatically with the operation of instruments like the Square-Kilometer- Array (SkA). The 
sensitivity of SkA will be around 100 times higher than the VLA sensitivity; it is expected that around 
20,000 new pulsars (including pulsars around black holes) will be discovered with SkA. The number of 
millisecond pulsars is expected to go up from its present value by a factor of 100. The initial operations 
of SkA will start around 2016, and the final operations are expected to begin around 2020. 

2 Bounds on the Nuclear Equation of State 

In 1939, Tolman, Oppenheimer and Volkoff performed the first neutron star calculations, assuming that 
such objects are entirely made of a gas of non- interacting relativistic neutrons [121 13] • The EoS of 
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Figure 3: Neutron star mass versus central energy density (eo = 140 MeV/fm 3 



such a gas is extremely soft, i.e. very little additional pressure is gained with increasing density, as can 
be seen from Fig. El and predicts a maximum neutron star mass of just 0.7 M (Figs. [3] and HI) at an 
unrealistically high density of 17 times the density of nuclear matter (Fig. [3]). It is interesting to note 
that the inclusion of interactions among the neutrons increases the star's maximum mass from 0.7 M 
to around 3 M & (Figs. [3] and H]). However, the radii of the latter stars are so big that mass shedding from 
the equator occurs for most stars of the sequence at rotational frequencies that are considerably smaller 
than 716 Hz observed for PSR J1748-2446ad, or B1937+21, 630 Hz (1.58 ms) [ID]. An interacting 
neutron gas thus fails to accommodate the rotational frequencies of the most rapidly rotating, observed 
neutron stars. The other extreme, a non-interacting relativistic neutron gas, fails too since it does 
not accommodate the Hulse- Taylor pulsar (M = 1.44 M ) [TJ], and also conflicts with the average 
neutron star mass of 1.350 ± 0.004 M derived by Thorsett and Chakrabarty [15] from observations of 
radio pulsar systems. More than that, recent observations indicate that neutron star masses may be 
as high as around 2 M Q . Examples of such very heavy neutron stars are Mj 0751+1807 = 2.1 ± 0.2 M 
US], M 4U1636+ 536 = 2.0 ±0.1 M [17], Mvei a x-i = 1.86 ± 0.16 Af Q [H], M CygX _ 2 = 1.78 ± 0.23 M 
[T9l [20] . Large masses have also been reported for the high-mass X-ray binary 4U 1700-37 and the 



2.44 ± 0.27 M [2T] and 



compact object in the low-mass X-ray binary 2S0921-630, M 4 u 1700-37 
M 2S0 92i-630 = 2.0 - 4.3M e 







], respectively. The latter two objects may be either massive neutron 
stars or low-mass black holes with masses slightly higher than the maximum possible neutron star mass 
of ~ 3M . This value follows from a general, theoretical estimate of the maximal possible mass of a 
stable neutron star [23J. If either one of the two objects 4U 1700-37 or 2S0921-630 were a black hole, it 
would confirm the prediction of the existence of low-mass black holes [23]. Conversely, if these objects 
were massive neutron stars, their high masses would severely constrain the EoS of dense nuclear matter. 



3 Nuclear Many-Body Models 



A vast number of models for the equation of state (EoS) of neutron star matter has been derived in the 
literature over the years. The majority of these models belong to either one of the following categories: 
1. Thomas- Fermi based models [251126] such as the new Thomas- Fermi approach of Myers and Swiatecki 
[25] . It is based on a Seyler-Blanchard potential generalized by the addition of one momentum dependent 
and one density dependent term 
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The upper (lower) sign in Eq. (JTJ corresponds to nucleons with equal (unequal) isospin. The quantities 
kp , T (= kp o /2m), and po are the Fermi momentum, the Fermi energy and the saturation density of 
symmetric nuclear matter. The potential's radial dependence is described by a Yukawa-type interaction 
of the form Y(r\2) = (47ra 3 ) -1 exp(— r^/ a) / {r^/a). Its strength depends both on the magnitude of the 
particles' relative momentum, P12, and on an average of the densities at the locations of the particles. 
The parameters £ and ( were introduced in order to achieve better agreement with asymmetric nuclear 
systems. The behavior of the optical potential is improved by the term <r(2p/po) 2/ ' 3 with the average 
density defined as p 2 / 3 = (p 2 / 3 + p 2 . )/2, where pi and p 2 the densities of interacting neutron or protons 
at points 1 and 2. The seven free parameters of the theory are adjusted to the properties of finite nuclei, 
the parameters of the mass formula, and the behavior of the optical potential [23] . 

2. Schroedinger-based models [31 [7J [271 1211 [29] are derived from Hamiltonians of the form 

^ = E^v 2 + E^+ E %, (2) 

i i<j i<j<k 

where V^- and Vijk denote two and three-nucleon interactions. The many-body equations are then 
solved, for instance, in the framework of the hole- line expansion (Brueckner theory), coupled cluster 
method, self-consistent Green functions technique, or variational approach. 

3. Relativistic nucler field-theoretical treatments such as relativistic mean field (RMF), Hartree-Fock 
(RHF), standard Brueckner-Hartree-Fock (RBHF), density dependent RBHF (DD-RBHF) [H [21 EDI ED 
1321 f33l 1341 135] which are based on a Lagrangian of the form C = Cb + Cm + C-int + Ci ept , where 

C B = ^(z 7 ^-m)V, (3) 
C M = W (d^d^ -m^)-^ £ (\F$F^ -mlAfA^) , (4) 

/ i=cr,8 A n=UJ,p \ Z ' 

Here, Cb and Cm are the free baryonic and the free mesonic Lagrangians, respectively, and interac- 
tions are described by Ci nt , where Fjj' = 9^4 — d u A^ is the field strength tensor of one of the 
vector mesons (k = u,p). In RMF, RHF and RBHF the meson-baryon vertices T a (a = a,u,5,p) are 
density-independent quantities which are given by expressions like r CT = ig a for the scalar o meson, 
= 9ul^ + (V2)(/^/2m)c?A[7 A , 7 M ] for uj mesons, etc. [2]. In the framework of the DD-RBHF scheme, 
the meson-baryon vertices T a are not only determined by Dirac matrices but depend on the baryon 
field operators ip [53] . 

4. The Nambu-Jona-Lasinio (NJL) model was originally introduced to describe nucleons with dynam- 
ically generated masses. In recent years it has become a popular model to describe quarks-hadron 
matter, and to explore the condensation patter of color superconducting quark matter J36[ [37J [381 SHI 
SOI SD S21 S3] • The flavor SU(2) NJL-Lagrangian, for instance, which has been applied in (43] to derive 
both nuclear matter and quark matter phases, is given by 

C = ^ - m q )tfj + G-K ((^) 2 - fcrV) 2 ) - GUVVY) 2 - G p (^t^) 2 

+G s {i, 1 ,Cr 2 (3 A ^ T )^ T C- 1 lb r 2 (3 A ^) . (6) 

Here m q is the current quark mass, ip is the flavor SU(2) quark field, and the coupling constants G n , 
G^ and G p characterize the qq interactions in the scalar, pseudoscalar and vector meson channels, while 
G s refers to the interaction in the scalar diquark channel [3SJ [43]. 

5. Aside from the NJL model, there are several other phenomenological models based on quark degrees 
of freedom, such as the quark meson coupling model, the cloudy bag model, the quark mean field 
model, and the chiral SU(3) quark mean field model [44]. In the latter model, quarks are confined 
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Figure 4: Mass-radius relationship of neu- 
tron stars and strange stars [6]. 



Figure 5: Mass-radius relationship of elec- 
trically charged strange stars. 



within baryons by an effective potential. The quark-meson interaction and meson self-interaction are 
based on SU(3) chiral symmetry. The chiral SU(3) quark mean field model was applied recently to the 
study of neutron stars and strange stars [44] . 

A collection of equations of state computed for several of these models is shown in Fig. [2j Mass-radius 
relationships of neutron stars based on these EoS are shown in Figs. [4] and [51 Strange star sequences 
are shown in these figures too. The strange star sequences shown in Fig. [5] show the impact of ultra- 
strong electric fields of the mass-radius relationship of strange stars. The energy-momentum tensor 
of such stars consists of the usual perfect-fluid term, which is supplemented with the electromagnetic 
energy-momentum tensor [45], TX = (p + pc 2 )u v u^ + pS^ + [F^ l F u i + 5 v fJ, FkiF kl /4 , K]/4:ir, where is the 
fluid's four- velocity, p and pc 2 = e are the pressure and energy density, respectively, and F^ v satisfies the 
covariant Maxwell equation, \{—g) 1 / 2 F u ^]^ = AnJ u (—g) 1 l 2 . The total mass of the star, contained within 
a radial distance r from the star's center, is given by dm(r)/dr = (4irer 2 )/c 2 + (Q(r)/c 2 r)(dQ(r)/dr). 
This relation shows that, in addition to the standard term originating from the EoS of the stellar fluid, 
the electric field energy too contributes to the star's total mass. The Tolman-Oppenheimer-Volkoff 
(TOV) equation of electrically charged stars is given by 

dp 2G [m(r) + *£(p- , + ^ Q(0 dQ{r) m 



dr ~ c 2 r 2 (l _ 2Gm(r) GQ 2 (r) \ ^ > 4^4 dr 

\ c 2 r r 2 c 4 / 



As shown in [45], electric fields can substantially alter the structure of compact stars. This is specifically 
the case for the masses of neutron stars, provided they posses net electric charges. For strange stars, 
however, even the maximum possible electric fields, ~ 10 18 V/cm, modifies the mass-radius relationship 
only very weakly as can be seen from Fig. [51 



4 Building Blocks of Neutron Star Matter 

4.1 Hyperons and baryon resonances 

At the densities in the interior of neutron stars, the neutron chemical potential, p n , is likely to exceed 
the masses, modified by interactions, of E, A and possibly 5 hyperons [2, [46]. Hence, in addition 
to nucleons, neutron star matter may be expected to contain significant populations of strangeness 
carrying hyperons. The thresholds of the lightest baryon resonances (A~, A , A + , A ++ ) are not reached 
in relativistic mean-field (RMF) calculations. This is different for many-body calculations performed 
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Figure 6: Composition of standard neu- Figure 7: Same as Fig. [61 but computed 
tron star matter, i.e. T = MeV and for the RBHF approximation |47j . 
Yl = for the RMF approximation. 



for the relativistic Brueckner-Hartree-Fock (RBHF) approximation where A's appear rather abundantly 
[17] . compare Figs. Eland CD Depending on the star's mass, the total hyperon population can be very 
large [16], which is illustrated graphically in Figs. [SJ and [9] for rotating neutron stars whose EoS is 
computed in the framework of the relativistic DD-RBHF formalism |33j . The stars shown in these 
figures have rotational frequencies ranging from zero to the mass shedding frequency, u^, which is the 
maximum frequency a star can have before mass loss at the equator sets in. This frequency sets an 
absolute upper limit on stable rapid rotation and follows from a metric of the form [21 SB] 

ds 2 = -e 2u dt 2 + e 2 ^ (d<f) - oo dtf + e 2 » d9 2 + e 2X dr 2 , (8) 

where u, ip, [i and A denote the metric functions which depend on the radial coordinate r, polar angle 
9 and, implicitly, on the star's angular velocity Q as well as on the angular velocity uj at which the 
local inertial frames are dragged along in the direction of the star's rotation. In Newtonian mechanics, 
mass shedding is determined by the equality between centrifuge and gravity, and is readily obtained as 
27rz/K = ^Jm/R 3 . Its general relativistic counterpart, computed here, is given by [21 US] 

The primes denote derivatives with respect to the Schwarzschild radial coordinate. Finally, we take 
a brief look at the composition of proto-neutron star matter. The composition of such matter is 
determined by the requirements of charge neutrality and equilibrium under the weak processes, Bi — > 
B 2 + / + Ui and B 2 + I —>■ Bi + z/;, where B\ and B 2 are baryons, and I is a lepton, either an electron 
or a muon. For standard neutron star matter, where the neutrinos have left the system, these two 
requirements imply that Q = J2i Qi^Bi +J2i= e ,/i Qi n i = (electric charge neutrality) and fi Bi = bi\i n — q^ 
(chemical equilibrium), where qi/i denotes the electric charge density of a given particle, and (n/) is 
the baryon (lepton) number density. The subscript i runs over all the baryons considered. The symbol 
H Bi refers to the chemical potential of baryon i, bi is the particle's baryon number, and qi is its charge. 
The chemical potential of the neutron is denoted by /i n . When the neutrinos are trapped, as it is the case 
for proto-neutron star matter, the chemical equilibrium condition is altered to \i Bi = bifi n — qi{jJ — ji Vl ) 
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Figure 8: Hyperon composition of a rotat- Figure 9: Same as Fig. [8] but in polar 
ing neutron star in equatorial direction. direction. 



and fi e — n Ve = /i M — where p, Vl is the chemical potential of the neutrino v\. In proto- neutron 
star matter, the electron lepton number Yl = (n e + n Ue )/riB is initially fixed at a value of around 
Y Le = Y e + Y Ue ~ 0.3 — 0.4 as suggested by gravitational collapse calculations of massive stars. Also, 
because no muons are present when neutrinos are trapped, the constraint Y Lfi = Y^ + Y UmU = can 
be imposed. Figures [10] and [11] show sample compositions of proto-neutron star matter and standard 
neutron star matter (no neutrinos) computed for the relativistic mean-field approximation. The presence 
of the A particle in (proto) neutron star matter at finite temperature is striking. As already mentioned 
at the beginning of this section, this particle is generally absent in cold neutron star matter treated in 
the relativistic mean-field approximation. 



4.2 Meson condensation 

The condensation of negatively charged mesons in neutron star matter is favored because such mesons 
would replace electrons with very high Fermi momenta. Early estimates predicted the onset of a nega- 
tively charged pion condensate at around 2uq. However, these estimates are very sensitive to the strength 
of the effective nucleon particle-hole repulsion in the isospin T = 1, spin S = 1 channel, described by 
the Landau Fermi-liquid parameter g', which tends to suppress the condensation mechanism. Measure- 
ments in nuclei tend to indicate that the repulsion is too strong to permit condensation in nuclear matter 
[4~9"1 [50] . In the mid 1980s, it was discovered that the in-medium properties of mesons may be 

such that this meson rather than the ir~ meson may condense in neutron star matter [5T| 152] . The 
condensation is initiated by the schematic reaction e~ — ► K~ +v e . If this reaction becomes possible in 
neutron star matter, it is energetically advantageous to replace the fermionic electrons with the bosonic 
K~ mesons. Whether or not this happens depends on the behavior of the K~ mass, m* K -, in neutron 
star matter. Experiments which shed light on the properties of the K~ in nuclear matter have been 
performed with the Kaon Spectrometer (KaoS) and the FOPI detector at the heavy-ion synchrotron 
SIS at GSI |53j. An analysis of the early K~ kinetic energy spectra extracted from Ni+Ni collisions 
showed that the attraction from nuclear matter would bring the K~ mass down to m* K - ~ 200 MeV at 
densities ~ 3uq. For neutron-rich matter, the relation m* K ~ /ttik- — 1 — 0.2n/rio was established [54], 
with rriK = 495 MeV the K~ vacuum mass. Values of around m* K - ~ 200 MeV may be reached by the 
electron chemical potential, /z e , in neutron star matter [21 HE] so that the threshold condition for the 
onset of K~ condensation, /i e = m* K , might be fulfilled for sufficiently dense neutron stars, provided 
other negatively charged particles (S~, A - , d and s quarks) are not populated first and prevent the 
electron chemical potential from increasing with density. We also note that K~ condensation allows 
the conversion reaction n — >• p + K~ . By this conversion the nucleons in the cores of neutron stars can 
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Figure 10: Composition of hot (T = Figure 11: Same as Fig. flOl but for stan- 
40 MeV) proto-neutron star matter for dard {Yl = 0) neutron star matter [57] . 
Y L = 0.3 [57]. 



become half neutrons and half protons, which lowers the energy per baryon of the matter. The relative 
isospin symmetric composition achieved in this way resembles the one of atomic nuclei, which are made 
up of roughly equal numbers of neutrons and protons. Neutron stars are therefore referred to, in this 
picture, as nucleon stars. The maximum mass of such stars has been calculated to be around 1.5 M & 
[55] . Consequently, the collapsing core of a supernova, e.g. 1987A, if heavier than this value, should 
go into a black hole rather than forming a neutron star, as pointed out by Brown et al. [21]. This 
would imply the existence of a large number of low-mass black holes in our galaxy [24] . Thielemann 
and Hashimoto [56J deduced from the total amount of ejected 56 Ni in supernova 1987A a neutron star 
mass range of 1.43 — 1.52 M . If the maximum neutron star mass should indeed be in this mass range 
(~ 1.5 M Q ), the existence of heavy neutron stars with masses around 2 M (Sect. [2]) would be ruled 
out. 



4.3 H-matter and exotic baryons 

A novel particle that could be of relevance for the composition of neutron star matter is the H-dibaryon 
(H=([w<i][<is][s'u])), a doubly strange six-quark composite with spin and isospin zero, and baryon number 
two [58J. Since its first prediction in 1977, the H-dibaryon has been the subject of many theoretical and 
experimental studies as a possible candidate for a strongly bound exotic state. In neutron star matter, 
which may contain a significant fraction of A hyperons, the A's could combine to form H-dibaryons, 
which could give way to the formation of H-dibaryon matter at densities somewhere above ~ 4n [59J. 
If formed in neutron stars, however, H-matter appears unstable against compression which could trigger 
the conversion of neutron stars into hypothetical strange stars [601 ED] ■ Another particle, referred to as 
an exotic baryon, of potential relevance for neutron stars, could be the pentaquark, Q + ([ud] 2 s), with 
a predicted mass of 1540 MeV. The pentaquark, which carries baryon number one, is a hypothetical 
subatomic particle consisting of a group of four quarks and one anti-quark (compared to three quarks 
in normal baryons and two in mesons), bound by the strong color-spin correlation force (attraction 
between quarks in the color 3 C channel) that drives color superconductivity [62]. The pentaquark 
decays according to 6 + (1540) — > K + [su] + n[udd\ and thus has the same quantum numbers as the 
K+n. 
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Figure 12: Sample composition of chemically equilibrated quark-hadron (hybrid star) matter 
function of baryon number density [21 [6] . 



4.4 Quark deconfinement 

It has been suggested already several decades ago [631 EH ESI [661 E3 EHl [691 EQ] that the nucleons in the 
cores of neutron stars may melt under the enormous pressures that exist in the cores, creating a new 
state of matter know as quark matter. From simple geometrical considerations it follows that quark 
confinement could occur at densities somewhere between around 2 — 10 n$. Depending on rotational 
frequency and neutron star mass, densities greater than two to three times no are easily reached in the 
cores of neutron stars so that the neutrons and protons in the cores of neutron stars may indeed be 
broken up into their quarks constituents pQ El El [71]. More than that, since the mass of the strange 
quark is only around 150 MeV, high-energetic up and down quarks will readily transform to strange 
quarks at about the same density at which up and down quark deconfinement sets in. Thus, if quark 
matter exists in the cores of neutron stars, it should be made of the three lightest quark flavors. The 
remaining three quark flavors (charm, top, bottom) are way to massive to be created in neutron stars. 
For instance, the creation of charm quark requires a density greater than 10 17 g/cm 3 , which around 
100 times greater than the density reached in neutron stars. A stability analysis of stars with a charm 
quark population reveals that such objects are unstable against radial oscillations and, thus, can not 
exist stably in the Universe [2J, [6] - The same is true for ultra-compact stars with unconfined populations 
of top and bottom quarks, since the pulsation eigen-equations are of Sturm-Liouville type. 

The phase transition from confined hadronic (H) matter to deconfined quark (Q) matter is char- 
acterized by the conservation of baryon charge and electric charge. The Gibbs condition for phase 
equilibrium then is that the two associated chemical potentials, /i n and /i e , and the pressure in the two 
phases be equal [71], Ph(//\ /i e , {ip x }, T) = Pq(//\ fi e , T), where Ph denotes the pressure of hadronic 
matter computed for a given hadronic Lagrangian C M ({ip x }), with {ip x } the field variables and Fermi 
momenta that characterize a solution to the field equations of confined hadronic matter, 

{i^d,-m x )^ x {x) = Fm x M(x)^ x (x), (10) 

M=<J,U1,7T,... 

R + m>(x) = ]T f^,(^(i), etc. (11) 

X=p,n,E,... 

The pressure of quark matter, Pq, is obtainable from the bag model. The quark chemical potentials 
fi u , fi d , /I s are related to the baryon and charge chemical potentials as /i u = (fi n — 2/i e )/3 and /i d = 
fi s = (fi n + /i e )/3. The Gibbs condition is to be supplemented with the global relations for conservation 
of baryon charge and electric charge within an unknown volume V containing A baryons. The first 
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r (km) r (km) r (km) 

Figure 13: Dependence of neutron star composition on spin frequency, u, for three sample 
compositions (left: hyperon composition, middle: quark- hybrid composition, right: quark- 
hybrid composition with quark matter in the color-flavor locked phase). The non-rotating 
stellar mass in each case is 1.4 M & . uk denotes the Kepler (mass-shedding) frequency. 



one is given by n = A/V = (1 — x) n n{u n , /A T) + xnq(/i n , /i e , T), where x = Kj/^ denotes the 
volume proportion of quark matter, Vq, in the unknown volume V, and and uq are the baryon 
number densities of hadronic matter and quark matter. Global neutrality of electric charge within 
the volume V can be written as = Q/V = (1 — x)?h /A T) + xgQ(/i n , T) + q^, with qt the 
electric charge densities of hadrons, quarks, and leptons. For a given temperature, T, these equations 
serve to determine the two independent chemical potentials and the volume V for a specified volume 
fraction x of the quark phase in equilibrium with the hadronic phase. After completion Vq is obtained 
as Vq = xV ■ The chemical potentials depend on the proportion x °f the phases in equilibrium, 
and hence so do all properties that depend on them, i.e. the energy densities, baryon and charge 
densities of each phase, and the common pressure. For the mixed phase, the volume proportion of 
quark matter varies from < % < 1 and the energy density is the linear combination of the two phases 
e = (1 — x)e H (/i r \ /i e , {V'x}' T) + X e Q,{^ n i /A T). Model neutron star compositions computed within the 
framework described just above are shown in Figs. [T2] and [T3J Possible astrophysical signals associated 
with quark deconfinement, the most striking of which being "backbending" of isolated pulsars, are 
discussed in [H [21 U [72] . 

4.5 Color-superconductivity of quark matter 

There has been much recent progress in our understanding of quark matter, culminating in the discovery 
that if quark matter exists it ought to be in a color superconducting state [251 E3, EH ES]- This is 
made possible by the strong interaction among the quarks which is very attractive in some channels 
(antisymmetric antitriplet channel). Pairs of quarks are thus expected to form Cooper pairs very 
readily. Since pairs of quarks cannot be color-neutral, the resulting condensate will break the local 
color symmetry and form what is called a color superconductor. The phase diagram of such matter is 
expected to be very complex [TSJ [74] . The complexity is caused by the fact that quarks come in three 
different colors, different flavors, and different masses. Moreover, bulk matter is neutral with respect to 
both electric and color charge, and is in chemical equilibrium under the weak interaction processes that 
turn one quark flavor into another. To illustrate the condensation pattern briefly, we note the following 
pairing ansatz for the quark condensate {tp^C^^il) 13 ^) ~ Aie^e/^i + A 2 e a/32 e/ a / b 2 + A3e a/33 e/ a / b 3, where 
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ipf a is a quark of color a = (r,g,b) and flavor f a = (u,d,s). The condensate is a Lorentz scalar, 
antisymmetric in Dirac indices, antisymmetric in color, and thus antisymmetric in flavor. The gap 
parameters Ai, A 2 and A 3 describe d-s, u-s and u-d quark Cooper pairs, respectively. The following 
pairing schemes have emerged: At asymptotic densities (m s — > or /i — > oo) the ground state of 
QCD with a vanishing strange quark mass is the color-flavor locked (CFL) phase (color-flavor locked 
quark pairing), in which all three quark flavors participate symmetrically. The gaps associated with 
this phase are A3 ~ A2 = Ai = A, and the quark condensates of the CFL phase are approximately 
of the form (ip^C^ip^) ~ Ae a/3X €f a f b x, with color and flavor indices all running from 1 to 3. Since 

e a ^ x ef a f b x = fif a $f b ~ fif b fif a one sees that the condensate involves Kronecker delta functions that link 
color and flavor indices. Hence the notion color-flavor locking. The CFL phase has been shown to be 
electrically neutral without any need for electrons for a significant range of chemical potentials and 
strange quark masses [76J. If the strange quark mass is heavy enough to be ignored, then up and 
down quarks may pair in the two-flavor superconducting (2SC) phase. Other possible condensation 
patterns are CFL-K , CFL-K + and CFL-7r°'~, gCFL (gapless CFL phase), ISC (single-flavor-pairing), 
CSL (color-spin locked phase), and the LOFF (crystalline pairing) phase, depending on m s , fi, and 
electric charge density (for references, see [B]). Calculations performed for massless up and down quarks 
and a very heavy strange quark mass (m s — > 00) agree that the quarks prefer to pair in the two-flavor 
superconducting (2SC) phase where A 3 > , and A 2 = Ai = 0. In this case the pairing ansatz 
reduces to (ipf a C , y 5 ipf b ) °t Ae a t,e a/33 . Here the resulting condensate picks a color direction (3 or blue 
in the example above), and creates a gap A at the Fermi surfaces of quarks with the other two out 
of three colors (red and green). The gapless CFL phase (gCFL) may prevail over the CFL and 2SC 
phases at intermediate values of rv? s j\i with gaps given obeying the relation A 3 > A 2 > A x > 0. 
For chemical potentials that are of astrophysical interest, fi < 1000 MeV, the gap is between 50 and 
100 MeV. The order of magnitude of this result agrees with calculations based on phenomenological 
effective interactions [751 EI] as we U as with perturbative calculations for \i > 10 GeV [TS]. We also 
note that superconductivity modifies the EoS at the order of (A//i) 2 [T9"l [50] , which is even for such 
large gaps only a few percent of the bulk energy. Such small effects may be safely neglected in present 
determinations of models for the EoS of quark-hybrid stars. There has been much recent work on how 
color superconductivity in neutron stars could affect their properties [73J [7U EH [821 [831 EI]. These 
studies reveal that possible signatures include the cooling by neutrino emission, the pattern of the 
arrival times of supernova neutrinos, the evolution of neutron star magnetic fields, rotational stellar 
instabilities, and glitches in rotation frequencies. 

4.6 Absolutely stable strange quark matter 

It is most intriguing that for strange quark matter made of more than a few hundred up, down, and 
strange quarks, the energy of strange quark matter may be well below the energy of nuclear matter, 
E/A — 930 MeV, which gives rise to new and novel classes of strange matter objects, ranging from 
strangelets at the low baryon-number end to strange stars at the high baryon number end [21 [6_1 [851 
I5B1 1ST] . The presence of electrons in strange quark matter may lead to the formation of an electric 
dipole layer on the surface of strange matter, with huge electric fields on the order of 10 19 V/cm. This 
peculiar feature enables strange quark stars to be enveloped in nuclear crusts made of ordinary atomic 
nuclei [851 EH EH]. Sequences of compact strange stars with and without (bare) nuclear crusts are 
shown in Fig. [H Since the nuclear crust is gravitationally bound to the quark matter, the mass-radius 
relationship of strange stars with crusts resembles the one of neutron stars and even that of white 
dwarfs [HD1 EI]- In contrast to neutron stars, however, strange stars obey M oc R 3 because they are 
self-bound and the mass density of quark matter is almost constant inside strange stars. The electrons 
surrounding strange quark matter are held to quark matter electrostatically. Since neither component, 
electrons nor quark matter, is held in place gravitationally, the Eddington limit to the luminosity that 
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a static surface may emit does not apply, and thus the object may have photon luminosities much 
greater than 10 38 erg/s. It was shown by Usov [92] that this value may be exceeded by many orders 
of magnitude by the luminosity of e + e~ pairs produced by the Coulomb barrier at the surface of a hot 
strange star. For a surface temperature of ~ 10 11 K, the luminosity in the outflowing pair plasma was 
calculated to be as high as ~ 3 x 10 51 erg/s. Such an effect may be a good observational signature of 
bare strange stars [92l |93| [9U [95]. If the strange star is enveloped by a nuclear crust however, which 
is gravitationally bound to the strange star, the surface made up of ordinary atomic matter would be 
subject to the Eddington limit. Hence the photon emissivity of such a strange star would be the same 
as for an ordinary neutron star. If quark matter at the stellar surface is in the CFL phase the process of 
e + e~ pair creation at the stellar quark matter surface may be turned off, since cold CFL quark matter 
is electrically neutral so that no electrons are required and none are admitted inside CFL quark matter 
[76] . This may be different for the early stages of a hot CFL quark star [96J. 



5 Latent Heat of Phase Transitions 

A neutron star is born with enormous reserves of rotational and thermal energy, which it looses over 
millions of years through processes such as magnetic dipole radiation from the magnetosphere, thermal 
radiation and electron winds from the surface and neutrinos from the core. In doing so the star cools 
down and changes structure as the radius diminishes with the rotation frequency-by as much as a 
few kilometers if it started out near its limiting Kepler frequency-shifting the interior boundaries of 
any phase transitions there in the process. Any latent heat evolved or absorbed by particles crossing 
these boundaries will contribute to the thermal evolution af the star, and in this section we discuss the 
formalism for describing these effects. The thermal evolution of a spherical star in General Relativity 
is given by a local energy balance equation [97J [98] : 

e-*-** d{L r e») _ _ ,Pdp _ x ^d(Te*) 

, 5 7\ — e nuc PC — — -+- e — , Li r — £ ±KKT C - , {LZ) 

4irr z or ot p at or 

where L r is the local luminosity given by the conduction of non-neutrino energy through a sphere of 
radius r, e nuc is the rate, per volume, at which nuclear reactions create non-nuclear energy, e„ is the 
neutrino emissivity per volume, p is the rest mass density, P is the pressure and II is the specific internal 
energy per unit rest mass, which in the absence of any phase transitions may be written 

LT = - / T c v (p,T')dT', (13) 
p Jo 

with cy the heat capacity per unit volume and T the temperature. Since P = —(dE/dV)T we then 
have dU/dp = P/p 2 so in this case 

dU dUdT dUdp c v OT P dp 

— = 1 = — 1 (1A) 

dt OT dt dp dt p dt p 2 dt ' K J 

which gives the usual heat balance equation for neutron stars 

d(L r e^) OT 

— e v — e cy— . (15) 



Anr 2 dr ' dt 

However in the presence of any phase transitions the heat capacity may be discontinuous and this will 
give rise to extra terms as particles cross the phase boundaries in either temperature for superfluid 
transitions or in density if there is a transition from hadronic matter to quark matter. To see this for 
a sharp phase transition with no mixed phase, we may write n as 



n = p" 1 



0(l-A/A Cil ) / T c v ,idr + e(A/A Ci2 -l) / r c v , 2 dT' 
Jo Jo 



(16) 
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where is the Heaviside step function, A can be either T, p (which may be discontinous at the phase 
transition) or some other convenient quantity depending on the nature of the phase transition, A c is the 
critical A for the phase transition and numerical subscripts identifies a variable with either of the two 
phases. Differentiating with respect to time we then get the same terms as above and additional terms 
from the Heaviside functions: 



~dt 



cydT P_dp 1 
P 



p dt p 2 dt 



6(X/X 



c,2 



l)dX 



A, 



c,2 



dt Jo 



cv, 2 dr' 



6(X/X cA - 1) OX 
dt 



X 



c,l 



CyidT' 



(17) 



Simplifying this we note that cy = T(ds/dT)y where s is the entropy density, integrate by parts and 
use the phase equilibrium conditions that P 1 = P2,Ti = T 2 , p± = p 2 to get 



dU _ c v or 
~dt ~ ~J~dt 



Pdp ldX^ 
- H T 

p 2 dt p dt 



<5(A/A C , 2 - 1) 5(A/A Cjl - 1) 



A 



-s 2 



c,2 



Ac, 



This is of course the well known result that each particle crossing a phase boundary contributes to the 
surrounding medium an energy [99] q = — T[a 2 — o"i] = —TAa, where a is the entropy per particle. If A 
is a continuous variable such as T or P the distinction between A Cj i and X C)2 is pointless and the delta 
functions may be taken outside, and a second order transition with continuous entropy density at the 
phase transition would then cause the last term in Eq. ffl8l) to vanish. 

A mixed phase may be treated similarly by replacing the Heaviside functions with a volume fraction 
X, such that 1 — x is the fraction of a given volume occupied by phase 1 and x is that occupied by phase 
2 so n = (1 — x)Ili + %n 2 and hence 



cM = cydT Pdp IdxdX r 

dt P dt p 2 dt P dx dt [S2 



(19) 



The transition to a superfluid state is a second order transition taking place through a mixed phase 
with Act oc —(1 — T/T c ) and lim-r^oX = 1 — (27rAo/T) 1//2 e _Ao//T for simple Fermi gasses [99]. The 
latent heat released in this process is customarily treated as an increase in the heat capacity letting 
cy — > cyR(T/T c ) with the effect of slightly slowing the cooling process - at temperatures around 0.2T C 
the reduction factor, R(T/T C ), goes to zero and has the opposite effect of significantly accelerating the 
cooling. This has been extensively treated in the literature and we refer to |100j for a recent review. 

The transition from hadronic to quark matter is a first order transition which will take place 
through a mixed phase [71] if finite size effects and charge screening do not prohibit such a phase 
|101[ 1102] (although see the comments in |103j ) . If the transition is sharp it happens at a specific 
critical pressure, where the density will be discontinuous and heat will be released there at a total rate 
of — e _</i TA(T(diV q /dt), with diVq/dt the number of particles making the transition per unit time. If the 
transition is through a mixed phase the pressure and average entropy density, s — (1 — x) s i + X s 2, vary 
continously in the region where the phase transition takes place, but the individual entropy densities will 
not be equal and particles making the transition will release heat at a local rate of —e~^(dx/dt)TAs, 
where x wm increase as the star becomes more dense. 

To estimate whether the contribution from these terms makes any significant difference or we are 
just carrying coal to Newcastle, we note that a transition from hadronic to quark matter may in a rough 
approximation be seen as a transition between free Fermi gasses at different densities-one of them with 
a bag constant in the energy density, but this does not affect the entropy density. Assuming both gasses 
are relativistic and degenerate the entropy per particle with Boltzmann's constant set to one is 



a 



7< 



(3tt 2 



3hc 



-Tn 



i 

" 3 ~ 



7c 0.02 



T 



MeV Vfm 



a 



(20) 
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where 7 C is the color degeneracy, n is the particle density which for quark matter is three times the 
baryon density, uqm, in the quark phase and in the hadronic phase is equal to the baryon density there, 
Tin- The resulting term in in Eq. f[T2l is then on the order of 



10 



T 2 

33^9 



fm~ : 



3n. 



QM 



fm 



dN q /dt _j 
-erg s 



10 57 /10 7 yr 



(21) 



where T 9 = T/10 9 K and dN q /dt is normalized to a whole star converted on a timescale of 10 7 years. For 
the quarks 7 C = 3 was canceled by a factor 1/3 since it takes 3 quarks to make a baryon. This is a rather 
small contribution unless the temperature is high or the star is changing structure fast. Furthermore the 
term is only positive if % < 3nQM, but this condition is of course subject to the specific assumptions 
we made here and would change in a more detailed treatment of the EoS. It should also be noted that 
this treatment accounts only for the latent heat from the deconfinement of the quarks, and that there 
would be further contributions from the subsequent weak reactions required to bring the matter into 
chemical equilibrium by converting <i-quarks to s-quarks. These reactions would contribute to the e nuc 
and e„ terms in Eq. (fT2"|) and must be treated similarly to the rotochemical heating discussed by [104J. 
As shown by |1U5] it is furthermore very important to know exactly where in the star a heating source 
is located, since sources located below the outer crust of the star tend to contribute only to the neutrino 
luminosity and make little difference to the surface temperature. In spite of such reservations one might 
still optimistically expect to find signals from the latent heat of a quark matter phase transition either 
early in a star's thermal history, when it is hot and spinning down rapidly, or from later episodes of rapid 
structure change. Such rapid structure change could result from glitches, which may be associated with 
the buildup and release of stress in a crystalline mixed quark hadron phase [106] . or the appearance of 
quark matter in a neutron star core, which softens the equations of state and allows the star to change 
structure rapidly possibly even through a core quake if a metastable state can be accessed |107] . Future 
work will explore these possibilities. 



6 Summary 

Obviously, our view of the interior composition of pulsars, which contain matter in one of the densest 
forms found in the Universe, has changed dramatically since their first discovery some 40 years ago. It 
has also become clear during that time period that all the ambient conditions that characterize pulsars 
tend to the extreme as well, rendering pulsars to almost ideal astrophysical laboratories for a broad 
range of physical studies. Owing to the unprecedented wealth of high-quality data on pulsars provided 
by radio telescopes, X-ray satellites-and soon the latest generation of gravitational-wave detectors-it 
seems within reach to decipher the inner workings of these enigmatic objects, and to explore the phase 
diagram of cold and ultra-dense hadronic matter from astrophysics. 
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